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• , ABSTRACT 

Context. BD-F20 1790 is a young active, metal-rich, late-type K5Ve star We have undertaken a study of stellar activity and kinematics 
for this star over the past few years. Previous results show a high level of stellar activity, with the presence of prominence-like 
^ , structures, spots on the surface and strong flare events, despite the moderate rotational velocity of the star In addition, radial velocity 

■ variations with a semi-amplitude of up to 1 km s"' were detected. 

^ ' Aims. We investigated the nature of these radial velocity variations, in order to determine whether they are due to stellar activity or 

the reflex motion of the star induced by a companion. 

Methods. We have analysed high-resolution echelle spectra, by measuring stellar activity indicators, and computing radial velocity 
(RV) and bisector velocity spans. Also two-band photometry was obtained to produce the light curve and determine the photometric 
^ ' period. 

CO \ Results. Based upon the analysis of the bisector velocity span, as well as spectroscopic indices of chromospheric indicators, like e.g. 

, Ca II H & K, Ha, and taking into account the photometric analysis, we report that the best explanation for the RV variation is the 

■ presence of a sub-stellar companion. The Keplerian fit of the RV data yields a solution for a close-in massive planet with an orbital 
^S) ' period of 7.78 days. The presence of the close-in massive planet could also be an interpretation for the high level of stellar activity 

detected. Since the RV data are not part of a planet search program, we can consider our results as a serendipitous evidence of a 
planetary companion. To date, this is the youngest main sequence star for which a planetary candidate has been reported. 
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1 . Introduction technique to detect exoplanets around young and active stars re- 

, . r , n , ,•• • quires, in addition, a careful characterization of stellar activity. 

H ; Smce the detection of the first planet orbiting a main sequence ^^^^^ -^^ ^^^jj^^. ^^^.^^^^ ^.^^^^^ 

.C^ star, 51 Peg (Mayor & Queloz 1995), the radial velocity (RV) ^j^^ ^ ^^^^^j j.^^^^^ ^ ,^,bsequent tem- 

method has become the most successful technique for detecting ^.^j ^^^^^^^^ RVs that may mimic a planetary reflex mo- 

exoplanets as the vast majority have thus far been discovered ^-^^ ^-^^ ^ -^^ j rotational period of the star (Saar 

in this way (Udry & Santos 2007). This method is especially ^ j^^^^^j^^^ ^^^^^ S^^^ ^^^^^ ^^j^^ planetary detections are 

efficient for giant planets in close-in orbits owing to the large ^ed by Queloz et al. (2001), Bouvier et al. (2007), Huerta 

radial velocities they indu ce in the host star. The use of the RV ^2008) and Huelamo et al. (2008). Thus the challenge in 

V, „ J r ~ 11 » J , ,u c ■ u using the RV technique to detect young planets lies in disentan- 

Based on observations collected at the German-Spanish 

Astronomical Center, Calar Alto, jointly operated by the Max- gl^g *e increased levels of stellar activity of young stars from 

Planck-Institut fur Astronomic Heidelberg and the Instituto de the RV signals of the planets. 

Astrofisica de Andalucfa (CSIC). Based on observations made with There is an absence of planets detected around stars younger 

the Italian Telescopio Nazionale Galileo (TNG) operated on the island j^an 100 Myr (Setiawan et al. 2007, Setiawan et al. 2008). Most 
of La Palma by the Fundacion Galileo Galilei of the INAF (Istituto searches for planetary companions have focussed mainly 

Nazionale di Astrofisica) at the Spanish Observatorio del Roque de . u i v . v. j f rt\r 

, , , r 1 T ■ J » ■ J ^ ■ V. J on stars older than 1 Gyr. Young stars were omitted from RV 

los Muchachos ot the Instituto de Astrolisica de Canarias. Based ., , i > m , ^ > 

on observations made with the Liverpool Telescope operated on surveys untd recently. Nevertheless, great effort has been made 

the island of La Palma by Liverpool John Moores University in the by several groups that have targeted young objects in their RV 

Spanish Observatorio del Roque de los Muchachos of the Instituto de searches of planetary companions. For example, surveys are 

Astrofisica de Canarias with financial support from the UK Science being carried out which focus on both nearby associations of 

and Technology Facilities Council. young stars and moving groups with ages ranging 10-500 Myr; 
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examples of which include p Pic (12 Myr), UMa association 
(300 Myr), Pleiades (100 Myr), IC 2391 (35 Myr), Hyades (700 
Myr), Taurus association (2 Myr), Chal (2 Myr), TWA (10 Myr) 
(Paulson et al. 2004, Paulson & Yelda 2007, Esposito et al. 2006, 
Huerta et. al 2007, Setiawan et al. 2007, Setiawan et al. 2008, 
Prato et al. 2008). Positive identification of planetary signatures 
from these eff'orts are few, with only two candidates to date: 
HD 70573 (Setiawan et al. 2007) and the controversial TW Hya 
(Setiawan et al. 2008). Planets orbiting around young stars are 
particularly valuable as they enable us to investigate some of the 
critical questions about the formation of both stellar and plane- 
tary systems. How, and at what stage planets form, what is the 
planet formation mechanism, and how they evolve are important 
questions which the study of young planetary systems will help 
to answer. In this paper we report strong evidence of a plane- 
tary candidate orbiting the young and active K5V star BD-i-20 
1790. Sect. 2, is an overview of the properties and our previous 
studies of this star. The observational strategy and data analysis 
are presented in Sect. 3. In Sect. 4 the nature of RV variations 
is investigated. An orbital solution for the data is presented in 
Sect. 5, and in Sect. 6 a discussion about planetary parameters, 
orbital solution, and how stellar activity and the planet are re- 
lated is shown. Finally, we summarize and offer some conclud- 
ing remarks in Sect. 7. 

2. BD+20 1790: An overview 

BD-H20 1790 was classified by Jeffries (1995) as a K5Ve star-, 
with a magnitude ofV- 9.9. Mason et al. (1995) identified this 
star as the optical counterpart for the 2RE J072343.6-H202500 
EUV source, located in the ROSAT All-Sky Survey. Lopez- 
Santiago et al. (2006) proposed its membership in the AB Dor 
kinematic moving group which has an estimated average age of 
50 Myr. By comparing the equivalent width of Li A 6708 A with 
the spectral type, Lopez-Santiago et al. (2006) derived an age 
estimate of 35-80 Myr. The main stellar parameters for BD-i-20 
1790 are compiled in Table[T] We obtained a value for the stellar 
radius from the measured rotational velocity and photometric pe- 
riod. Our estimated radius agrees with the previous K5 V spectral 
classification (from Carrol & Ostlie (2007) tables). Adopting this 
spectral type, we used the K5V temperature from the CaiTol & 
Ostlie (2007) tables. In conjunction with the photometric param- 
eters, this enabled us to derive the luminosity, mass and surface 
gravity. Errors in the parameters were estimated by following 
the method of propagation of errors, i.e., the uncertainties were 
calculated from the errors in the variables involved in the deter- 
mination of each parameter. It has been assumed null correlation 
between the different variables, in principle independent of each 
other. In order to test if assuming a fixed value for Teff has a 
non-negligible effect in the errors computation, we investigated 
whether an error in Teff could translate into uncertainties of de- 
rived parameters. We have considered an input error in reff~10 
K and conducted an analysis in the propagation of T^^ error. 
Based on this analysis, we have noticed that not consider the er- 
ror in the T^s leads to an underestimation in the errors in mass 
and log g, providing unreliable error bars for these parameters. 

The X-ray luminosity was calculated using the count rates 
and HRl hardness ratios from the ROSAT All-Sky Survey. By 
combining the conversion factor Ci, computed by the formula 
from Fleming (1995), and the distance estimated by Reid et al. 
(2004), the stellar X-ray luminosity was calculated as Lx - 1.6+ 
0.5 10^'^ ergs-\ 

We compute a preliminary value of metallicity by using a grid of 
Kurucz et al. (1993) ATLAS9 atmospheres and the 2002 version 



of MOOCQ synthesis code (Sneden 1973). Atmospheric mod- 
els were constructed with the data given in Table [1] We used 12 
Fe I lines selected from Gonzalez et al. (2001). We also calcu- 
lated a value of metallicity by using 7 Fe I lines in the MOOG 
Abfind routine. We find an average value of A[Fe] = 7.82+0.20 
which, when assuming a solar value of A[Fe\ - 7.52, results in a 
[FejH] - 0.30+0.20. As mentioned, this is a preliminary value, 
although compared with the average metallicity of stars of so- 
lar neighbourhood, we still could consider the star as metal-rich 
within the error bars. 

In a recent paper Carpenter et al. (2008) derived the tempera- 
ture, gravity and metallicity for BD+20 1790, being their values 
of Teff = 4408 K and log g - 4.50, very close to the correspond- 
ing values presented in Table [T] We also pointed out that the 
difference between metallicity values may be explained by the 
fact that Carpenter et al. (2008) only assumed a fixed metallicity 
of [FejH] - 0.0, but not actually compute it. 

The Li I abundances analysis was also done in standard local 
thermodynamic equilibrium (ETL) using MOOG and ATLAS9 
in the same way as with the metallicities. Abundances were 
derived by fitting synthetic spectra to the data. To determine 
Li abundances we perform a spectral synthesis around the Li 
1 6707 A resonance doublet, fitting all spectra between 6702 and 
6712 A, taking into account the relation between Li6 and Li7 
isotopes. We determine an average value of lithium abundance 
of logNiLi) = 1.03+0.04 (where logN{Li) = log (Li/H)+12). 

In order to study the stellar activity and the kinematics, we 
have carried out both spectroscopic and photometric monitoring 
over the past few years: high temporal and spectroscopic resolu- 
tion and two band photometry. The simultaneous study of pho- 
tospheric and chromospheric active regions is a powerful tool 
that allow us to trace, reconstruct and model the puzzle of the 
magnetic field topology, since these active regions are the finger- 
prints of magnetic fields (Collier Cameron 2001, Catalano et al. 
2002, Frasca et al. 2005, Collier Cameron et al. 2002). Strong 
chromospheric activity was detected in several observing runs, 
described by Hernan-Obispo et al. (2005, 2007). In spite of the 
fact that the rotational velocity is not very high, vsin/ ~ 10 km s"' 
(Lopez-Santiago et al. 2006), all activity indicators are in emis- 
sion above continuum, from Ca ii H & K, to Ca n IRT lines (see 

Fig. IB- 
Through the study of profile line asymmetries of Ha and Hj6 
lines, prominence-like structures have been detected in the chro- 
mosphere of the star (Hernan-Obispo 2005, 2007). These can be 
observationally detected as transient absorption features super- 
imposed on the line profile that are interpreted as the presence of 
cool material embedded in the surrounding hotter corona and co- 
rotating with star (Collier Cameron & Robinson 1989a,b, Collier 
Cameron & Woods 1992, Jeffries et al. 1993, Byrne et al. 1996, 
Elbe et al. 1998, Barnes et al. 2000, Donati et al. 2001). Several 
completed prominence-like transients have been detected with 
durations of orders of a few hours (see Hernan-Obispo 2005 
for details). Modeling these chromospheric phenomenae is an 
important challenge in this case, due to the detection of these 
prominence-like structures in unstable positions, far from equa- 
torial regions (Ferreira 2000, Jardine et al. 2001, Jardine & van 
Balegooijen 2006). 

In addition, strong large optical flare events were observed. 
The gradual decay of the flares was observed for up to 5 hours. 
Fig. [T] compares the activity indicators for the quiescent state 
and flare state. The energy released is on the order of ~10^^ erg. 



' The source code of MOOG 2002 can be downloaded at 
http ://verdi . as .utexas .edu/moog .html | 
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Table 1. Stellar Parameters of BD+20 1790 



Parameter 


Value 


Spectral Type 


K5 V 


B-V 


1.15 




r\ n'i 1 r\ r\r\ ji 

u.oJ + u.uy Mq 


^eff 


4410 K 


logg' 


4.53 ± 0.17 




110 ± 3 mA 


Distance'^ 


25.4 ±Apc 


Age- 


35 - 80 Myr 


vsin;'' 


1 A Al 1 A An 1^*^ .. — 1 

1U.U3 + 0.4/ km s 


pa 
phot 


2.801 ± 0.001 days 


50.41 degrees 




0.71 ±0.03 Rq 


[FelHf 


0.30 ± 0.20 


logN{Lif 


1.03 ±0.04 




1.6 ±0.5 10^' ergs-' 




0.17 ± 0.04 Lq erg S-' 



" This paper 

'° From Carrol & Ostlie, 2007 

- From Lopez Santiago et al. 2006 

From Lopez Santiago 2005 
"= From Reid et al. 2004 



while for largest solar flares the released energy is about ~10^^- 
10^^ erg, thus ranging the flares of BD+20 1790 on the so-called 
superflare regime (Rubenstein & Schaefer 2000). 

The photometric observations yielded a light curve with ev- 
idence of rotational modulation, the semi-amplitude of which is 
up to Ay~ 0."'06 and indicates the presence of spots on the sur- 
face. The period analysis of the entire set of observations reveals 
a photometric period of 2.801 (+ 0.001) days, in agreement with 
the period given by the SuperWASP photometric survey (Norton 
et al. 2007). 

A detailed and completed study of the chromospheric and 
photospheric activity characterization will be published in a 
forthcoming paper (Hernan-Obispo et al. 2009b, in prep.). 



3. Observations and Data Analysis 

In order to study and characterize active regions at photospheric 
and chromospheric levels, we carried out photometric and spec- 
troscopic observations of the target. 



3. 1 . Spectroscopic data 

The observational strategy was designed to spectroscopically 
monitor chromospheric activity indicators with high temporal 
and spectral resolution. High resolution echelle spectra were ob- 
tained during four observing runs, from 2004 to 2007, detailed 
in Table |2] The exposure times ranged from 900 s to 1200 s, 
depending on weather conditions, in order to obtain a S/N typ- 
ically greater than 140 for SARG runs and 80 for FOCES runs. 
The spectra in the time series observations were separated only 
by the CCD readout time, thus enabling us to obtain the high- 
est temporal resolution possible. Our initial temporal cadence 
was designed to detect prominence-like transient features in the 
Balmer lines. Spectral types and RV standards were acquired 
with the same setup and configuration as the target. These stan- 
dards were reduced and analysed in the same way as the tar- 




Wavelength (A) 




3934 
Wavelength (A) 




4862 
Wavelength (A) 



Fig. 1. Chromospheric activity indicators. The dashed line indi- 
cates quiescent state, while solid line indicates flare state. From 
top to bottom and left to right: He I D3 region, Ca 11 K, Ha and 



get. The data were bias-subtracted, overscan-corrected and flat- 
fielded using standard routines in IRAfQ package. 



' IRAF is distributed by the National Optical Observatory, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation. 
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The wavelength caHbration was obtained by taking spectra 
of a Th-Ar lamp. Using Coude spectrographs allowed a stable 
environment for the wavelength calibration, since flexures are 
not possible. Details about the spectrographs used can be seen 
in Pfeiff'er et al. (1998) for FOCES spectrograph and Gratton 
et al. (2001) for SARG spectrograph. In order to enhance the 
accuracy in calibration we used about 10-12 lines identified per 
order; across all orders for SARG spectra and about 80 orders 
for FOCES spectra. The orders were calibrated simultaneously 
and the total fit has an rms value typically lower than 0.003 A. 
The spectra were normalized by a polynomial fit to the observed 
continuum. 

Heliocentric radial velocities were determined using a 
weighted cross-correlation method. The spectra of the star were 
correlated order by order against spectra of several RV standards 
with similar spectral type. Orders with chromospheric features 
and telluric lines were excluded. We calculated the uncertainties 
based on the cross-correlation peak height and the antisymmetric 
noise as described by Tonry & Davis (1979). Also, by measur- 
ing RVs of the standard stars, we estimated the systematic errors 
and the accuracy of the RV measurements with our instrumental 
setup. The accuracy between standards for the same run and be- 
tween runs is less than 0.05 km s"'. 

Additional echelle data were acquired in DDT mode at the 
FOCES spectrograph on December 2008. The telescope config- 
uration and the setup were identically to previous FOCES runs, 
except for two nights in which a different CCD was used. Data 
were taken over 10 consecutive nights but due to bad weather 
conditions only five nights were acquired. Because of the time 
limitation in DDT mode, only one RV standard was observed. 

3.2. Photometric data 

The purposes of these observations were to determine the photo- 
metric period and to look for photometric variability. In addition, 
the study of the light curve, as well as spectrocopy, allow us to 
characterize the active regions in the photosphere (Catalano et 
al. 2002, Frasca et al. 2005, Biazzo et al. 2007). 

CCD differential aperture photometry was obtained using the 
2.0 m fully robotic Liverpool Telescope (Steele et al. 2004) at the 
Observatorio del Roque de los muchachos in La Palma, Spain. 
The observations were scheduled in monitoring mode. We ob- 
tained 22 photometric epochs during November and December 
2007. Our observational strategy permitted us to obtain a photo- 
metric epoch every 3 nights on average. Each epoch consisted 
in alternating r' and g' exposure^ thereby obtaining quasi- 
simultaneous two band photometry. Custom made softwarqjwas 
used to automatically extract the photometry. By analysing intra 
night scatter we can infer a photometric accuracy of 3 mmag and 
4 mmag per exposure (r' and g' bands respectively, see Fig. 
We fit the best sine-wave model to the photometry sampling 
many periods between 0.1 and 50 days, on both bands. Plotting 
the post-fit residuals as a function of the period, a very strong 
minimum on the post-fit residuals is found at 2.801 + 0.001 days 
in both bands (see Fig. |3]l. We note that the period and the am- 
plitude are similar with those given by the SuperWASP survey 
(Norton et al. 2007). 

The difi'erent amplitude in each band is consistent with large 
spot or spot group covering at least 4% of the surface. As can be 
seen in Fig. |2] the amplitude is larger at shorter wavelength, i. 



^ Sloan r' and g' filters were used 

" ATP, Automatic TATOOINE 

|http://www.am.ub.es/~anglada/atp/atp-testing.htmi 
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0.2 0.4 0.6 0.8 1 
Phase 

Fig. 2. Photometry phased to the 2.801 days period. A linear 
trend and a zero point have been subtracted to both bands The 
residuals with respect to a simple sine-wave model are show in 
the lower panel. 
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Fig. 3. Postfit residuals to the photometry as a function of the 
period. The sharper minima correspond to the 2.801 days period 
in both bands. The RVs period is marked in grey to illustrate the 
absence of related photometric signals. 



Photometry. 



e. at g' band in this case. This color variation is correlated with 
variation in magnitude. The star appears redder when fainter, at 
minimum light and therefore bluer when brighter, at maximum 
light. 

The full analysis of the photometry and its relation to the star 
activity requires simultaneous discussion with the spectroscopic 
data and a more detailed study of the star will be presented else- 
where (Hernan-Obispo et al. 2009b, in prep.) 
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Table 2. Observing runs 



Date 


Telescope 


Instrument 


CCD chip 


Spect. range 


Orders 


Dispersion 


FWHM'^^ 


N. Obs. 








# 


(A) 




(A/pix) 


(A) 






29/03-6/04 2004 


2.2m" 


FOCES 


2048x2048 24yum Site#ld 


3720 - 10850 


100 


0.04-0.13 


0.08- 


0.35 


19 


21-22/11/2004 


TNG'' 


SARG 


2048x4096 13.5yumEEV 


4620 - 7920 


52 


0.07-0.11 


0.07- 


0.17 


43 


15/04/2006 


TNG'' 


SARG 


2048x4096 13.5yumEEV 


4620 - 7920 


52 


0.07-0.11 


0.07- 


0.17 


14 


2-5/10/2007 


2.2m" 


FOCES 


2048x2048 24y[(m Site#ld 


3720 - 10850 


100 


0.04-0.13 


0.08- 


0.35 


10 


12-13/12/2008 


2.2m" 


FOCES 


2048x2048 15y[(m L0RAL#1 li 


3830 - 10850 


96 


0.03 - 0.07 


0.09- 


0.26 


2 


19-21/12/2008 


2.2m" 


FOCES 


2048x2048 24y[im Site#ld 


3620 - 7360 


100 


0.04-0.13 


0.08- 


0.35 


3 



" 2.2 m telescope at the German Spanish Astronomical Observatory (CAHA) (Almeria, Spain). 
'' 3.58 m Telescopio Nazionale Galileo (TNG) at Observatorio del Roque de los Muchachos (La Palma, Spain). 
The spectral resolution is determined as the FWHM at the arc comparison lines ranges. 



4. On the nature of the RV variations 

Variations in the RV peak-to-peak amplitude of up to ~2 km s"' 
were observed during all the observing runs. These variations 
are significantly larger than the individual measurement errors 
(0.10 to 0.20 km s"') or the systematic error (0.05 km s"'), even 
when we consider the scatter between runs with different spec- 
trographs and setups. 

4.1. Searching for periodical signals on RV 

A Least squares periodogram (see Appendix A) reveals one very 
significant peak at 7.783 days (see Fig.|4^). The data set contains 
91 independent RV measurements. However, many of them are 
clustered together within groups of a few hours. The values we 
use to generate the periodogram and for orbital fitting (shown in 
Table[3]l, are averaged on a nightly basis. Fig. |4}5 shows the em- 
pirical False Alarm Probability (FAP) as a function of the Power 
The 7.783 days peak has a FAP of 0.35%. 

It is worth noting that the RV period is larger than the pho- 
tometric period. Nevertheless to test if the RV period could arise 
from rotational modulation we searched for significant frequen- 
cies in the data points of the photometry. There is no significant 
power at the RV period, and no secondary peaks are found in 
the aliasing frequencies of the RVs or the photometric period 
after the main signals are removed. To illustrate the absence of 
related photometric signals, we marked the RV period in Fig.|3] 
that shows the post-fit residuals of photometric data. In addition 
to this, there is no signal at photometric period in the RV data, 
as can be seen in Fig.|4^, that shows the RV periodogram. 

4.2. Stellar activity jitter 

It is well known that spurious RV variations can be induced by 
stellar activity, especially due to changes in the profile of spectral 
lines caused by the presence of active regions, the so-called stel- 
lar jitter (Saar & Donahue 1997, Saar 2009). The high level of 
activity detected in BD4-20 1790, induced us at first to relate RV 
variations with active regions. Since we ruled out the possibility 
of variations due to systematic errors or any seasonal effect, the 
main concern was to determine if stellar activity was responsi- 
ble. 

It is widely accepted that the relationship of bisectors of the 
cross-correlation function (CCF) and RV is a powerful method 
to determine whether the RV variation may be due to stellar ac- 
tivity or a planetary companion (Queloz et al. 2001, Martmez- 
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Fig. 4. a. Up: Least-Squares Periodogram of the nightly av- 
eraged radial velocity measurements. The 7.78 days peak has 
a FAP of 0.35%. The dotted horizontal line illustrates a FAP 
lower than 1 % and the dashed horizontal line a FAP lower than 
5%.b. Down: Empirical False Alarm Probability as a function of 
the Power (red line). The gray bars illustrate the distribution of 
False Alarms with an arbitrary normalization used to derive the 
Empirical FAP. Note that the Y axis is in logarithmic scale. 



Fiorenzano et al. 2005). The CCF was determined by using the 
same procedure as for the RV case, computing it for the regions 
which include the photospheric lines which are more sensitive to 
spot presence, while excluding chromospheric lines and telluric 
lines. The bisector inverse slope (BIS), defined as the difference 
of the average values of the top and the bottom zones, was com- 
puted to quantify the changes in the CCF bisector shape by us- 
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Table 3. Radial Velocity 



JD days 


RV (km/s) 


tr (km/s) 


2452388.3341" 


9.23 


0.19 


2452389.3513" 


8.94 


0.14 


2452390.3670" 


8.52 


0.38 


2453099.3573" 


7.82 


0.06 


2453100.3692 


6.96 


0.10 


2453101.3748 


7.34 


0.07 


2453102.3876 


7.84 


0.05 


2454375.6480 


7.96 


0.08 


2454378.6804 


7.72 


0.04 


2453331.6400 


8.71 


0.03 


2453332.6800 


8.16 


0.03 


2453841.4250 


7.73 


0.03 


2454812.7429 


7.76 


0.21 


2454813.7240 


7.67 


0.18 


2454820.5057 


7.73 


0.16 


2454821.5126 


7.53 


0.16 


2454822.5483 


7.96 


0.14 



From Lopez-Santiago 2005 




7 7.5 a S.5 

RV (km/s) 



Fig. 5. Bisector velocity span vs. Radial velocity for all the 
observing runs. Symbols represent the different runs: stars for 
FOCES 04A, diamonds for SARG 04B, circles for SARG 06A 
and triangles for FOCES 07B. The lack of correlation indicates 
that RV variations are not due to stellar activity. 

ing the method described by Queloz et al. (2001). In choosing 
the span zones we avoided wings and cores of the CCF profiles, 
where errors of bisectors measurements are large. In Fig. |5] it 
can be seen that there is a lack of coiTelation between the BIS 
and RV variation for all the observing runs. This indicates that 
the RV variations are not due to variations in the asymmetry of 
the photospheric lines profile, and subsequently not due to stel- 
lar activity variations. The least squares periodogram of bisec- 
tors shows two tentative peaks around 2.8 days, and is shown in 

We have estimated the stellar jitter from Santos et al. (2000), 
that takes into account the Ca ii H & K index. Assuming an av- 
erage value for Ca n H & K index of about -4.2 and by using 
eq. [4], we derived a value for the stellar jitter of up to 10 m s"' . 
This stellar jitter is added in quadrature to the RV error. 
As an additional test we investigated the variation of stellar ac- 
tivity indicators, especially those that are ascribed to the pres- 
ence of plage-like structures on the chromosphere, like B aimer 
lines, Ca n H & K and Ca ii IRT. The emission flux for these 
lines in active stars usually shows a periodic modulation (and 




Period (days) 

Fig. 6. Periodogram for bisectors of all runs. There is no clear 
period for bisector variations. 



subsequently the spectroscopic indices) which is most likely due 
to rotational modulation of plage-like structure emission. As is 
shown in Sect. 2, all chromospheric activity indicators are in 
emission above the continuum, indicating a very high level of 
activity. To avoid the photospheric contribution to the spectral 
profiles we applied the spectral subtraction technique described 
in detail by Montes et al. (1995). This technique makes use 
of the program STARMOD developed at Penn State University 
(Barden 1985) and lately modificated by Montes et al. (1995). 
Also, in order to control the error and minimize the uncertain- 
ties, some routines of the astronomical data reduction package 
REiDmE developed at Universidad Complutense de Madrid 
(Cardiel 1999) were used. In these subtracted spectra, spectro- 
scopic indices have been defined and computed following Saar 
& Fisher (2000), Kiister et al.(2003), Bonfils et al. (2007). Both, 
Ca II IRT and Ca ii H & K indices were only determined for 
FOCES runs, due to the wavelength range coverage of the spec- 
trograph. To avoid contamination from telluric lines we only 
consider the 8662A Ca n IRT line. We searched for periodic 
signals in the spectroscopic indices by computing their Least 
squares periodograms. Fig. |7] shows the variation with time (or- 
bital phase folded in this case) for Ca ii IRT, Ca ii H & K, Ha and 
Hy6 indices. The corresponding periodogram computed shows 
more noise rather than a clear signal. This result is also seen 
in the indices figures as a non-modulation of the activity index. 
As an example, Fig.[8]shows the periodogram for the Ho- index. 

As pointed out by Walter (1994), the rotational modulation 
of chromospheric lines due to plages is not always detectable 
in very active stars. Furthermore, in this case the flares could 
contaminate the data, masking the actual period of variation of 
the indices. In order to investigate this possibility we removed 
the data affected by flare events. Due to the different wavelength 
range coverage of spectrographs, we considered only Ha and 
H/3 indices. For Ha index, we have found a tentative rotational 
modulation with a period of 2.77 days, similar to photometric 
period (see Fig.|9]l. However, the postfit residuals show in Fig.fTOl 
that this could be a misleading signal, even pure noise. For H/3 
index, no clear modulation has been found. 
The lack of variability of BIS and spectroscopic indices with RV 

^ http://www.ucm.es/info/Astrof/software/reduceme/reduceme.html 
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Fig. 7. Spectroscopic index for chromospheric activity indica- 
tors, phased folded orbital period. From top to bottom: Ha 
(squares), HyS (circles), Ca n IRT (triangles) and Ca ii H & K 
(stars). The dashed line is indicating the quiescent state. Error 
bars for indices are of order of 0.001 
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Fig. 8. Periodogram for Ha index. There is no clear period for 
index variations. 



CO 
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Fig. 9. Ha index for the data without flare events. It can be seen 
a modulation with a period of about 2.77 days, similar to photo- 
metric period. 




Period (days) 

Fig. 10. Postfit residuals to Ha index for no-flare data as a func- 
tion of the period. There is no clear period for index variations. 



4.3. RV wavelength dependence 

Desort et al. (2007) (hereafter D07) pointed out that the color de- 
pendence (with wavelength) of the RV peak-to-peak amplitude 
period, and the absence of a photometric period larger than 2.8 with spots can be used as a diagnostic to discriminate between 
days, strongly support the planetary companion hypothesis. stellar activity or planetary companions. Due to the contrast be- 
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tween spots and the surrounding photosphere being greater in 
the visible than at IR wavelengths, it is expected that an attenu- 
ation of RV amplitude towards red wavelengths would be seen. 
Observationally this effect has been shown by e.g. Martin et al. 
(2006), Huelamo et al. (2008) and Prato et al. (2008). If the 
RV variations are due to a planet, the RV amplitude should be 
the same at every wavelength range. We investigated a possi- 
ble chromatic dependence by computing the RV in two differ- 
ent ranges of wavelength, one for red and near-IR wavelengths 
(7650 to 10000 A) and the other for blue (4300 to 4800 A). The 
resulting RV peak-to-peak amplitude is 2.19+0.20 km s"' for 
the near-IR range and 2.20+0.20 km s ' for the blue range. The 
values differ by only 0.5% and agree within the uncertainties. 
Additional RV infrared follow-up can allow us to confirm this. 
In a forthcoming paper (Hernan-Obispo et al. 2009d, in prep.), 
we will present the first results of the study of the RVs of BD+20 
1790 in the near-IR range. 

4.4. RV variation by empirical spots and plages? 

To estimate an order of magnitude of the expected RV ampli- 
tude due to spots, we use empirical relations derived by Saar & 
Donahue 1997 (hereafter SD97) and D07. These relations con- 
nect the RV amplitude with the spot filling factor /j and vsini. 
We consider both relations by D07 and SD97, because D07 re- 
lations take into account the spectral type and the whole spectral 
range (except telluric and chromospheric lines) to compute the 
empirical RV, whereas SD97 uses a single line and G5V spectral 
type. Using SD97 eq. [1], we derived an amplitude of up to 575 
m s"' and by using D07 eq [5] we similarly estimate an ampli- 
tude of up to 600 m s"'. As mentioned, these results are taken 
as a quantitative estimation. There are more effects that are not 
taken into account here, like the spot location at stellar surface 
given by the colatitude 6, and the spot temperature. SD97 eq. 
[1] and D07 eq. [5] considered the simple case of an equatorial 
spot, but SD97 assumed a Tspoi = K and D07 assumed an spot 
temperature 1000 K cooler than the photosphere. The difference 
between the RV amplitude derived from both equations could be 
due to this different spot temperature. 

On the other hand, we can estimate the spot filling factor that 
could produce the RV signal of our data. We considered an aver- 
age semi-amplitude of 1 km s"'. The fs estimated from SD97 is 
therefore 23% while D07 indicates 19%. The fs measured from 
photometric variation is about 4%. These results indicate that the 
spot filling factor needed to explain the RV variation due purely 
to spots is not in agreement with the photometry. 

Saar (2003) and Saar (2009) showed significant efforts to 
model plage-induced RV jitter Although the models are mostly 
applicable to solar-like stars, we could estimate the plage fill- 
ing factor fp that could produce the RV signal by using the Saar 
(2009) equation that connects the RV amplitude with vsini > 6 
km s"^ This fp estimated is about 70%, that strongly suggets 
that the RV variation is not due to chromospheric plages. 

4.5. What would the RV signal be without a planet? 

It is important to remark that empirical relations derived by 
SD07 and D07 do not take into account the chromatic effect 
of spots on the RV signal. We therefore investigated how much 
RV signal would be expected in the absence of a planet, and the 
degree of RV attenuation with wavelength (assuming the RVs 
are due to cool spots). In order to quantify the attenuation if the 
cause of variations were spots, we try to investigate how much 




4000 6000 8000 10* 

Wavelength (A) 



Fig. 11. Radial velocity amplitude variation with wavelength, 
computed for LO Peg profiles 

spots affect the line profiles. However, BD+20 1790 has a low 
vsini to model the photosphere by generating Doppler imaging 
spot maps. To carry out a realistic approximation to the prob- 
lem, we construct realistic spot maps by using spectra of another 
star with similar characteristics, LO Peg, that is widely stud- 
ied in the literature and its photospheric activity is well-known 
(Jeffries & Jewell 1993, Jeffries et al. 1994, Elbe et al. 1998, 
Elbe et al. 1999, Barnes et al. 2005). LO Peg is a K5V-K7Ve 
star, identified by Jeffries & Jewell (1993) as a member of the 
Local Association, with an estimated age of 20-30 Myr Jeffries 
et al. (1994) determined the inclination to be 50°. The level of 
activity is similar to BD+20 1790, but LO Peg is a rapid rota- 
tor (vsim -69 km s"'). The LO Peg photometry suggests a spot 
filling factor of up to 1.5%. 

Using the Doppler imaging program, DoTS (Collier 
Cameron 1997), and an input starspot image derived for LO Peg 
(Barnes et al. 2005), we generated a set of line profiles for a 
star with vsini = 10 km/s (i.e. matching that of BD+20 1790) 
over a complete rotation phase. The profiles thus contain asym- 
metries due to starspots from the observed LO Peg image. We 
used appropriate temperatures for the BD+20 1790 photosphere 
and estimated the spots to possess temperatures which were up 
to 1000 K cooler Profiles were generated for the three differ- 
ent wavelengths of 4000 A, 6717 A and 10000 A. The radial 
velocity variations were then calculated in order to estimate the 
relative amplitudes due to spot induced variations at each of the 
three wavelengths. 

The RV attenuation with wavelengths relative to 4000 A is 
16% at 6717 A and about 30% at 10000 A , as illustrated in 
Fig-El Assuming at a first approach the same /, for LO Peg and 
BD+20 1790, the RV signal for BD+20 1790 should be about 1 .5 
kms-i at 10000 A. 

However, in Hernan-Obispo et al. 2009d (in prep.), we find 
only about 0.5% attenuation in the near-IR region relative to the 
visible, 6717 A region. This result is an additional argument in 
support the existence of the planetary companion. 

4.6. RV jitter from flares 

We have estimated the rate of flare occurrence as the fractional 
amount of the total observing time (for all runs) where a flare 
was detected. Thus, we get a flare frequency of occurrence of ~ 
40%. This higher rate raises the question of how much RV jitter 
we should expect from large flares, if any. Saar (2009) presents 
the first approach to this issue, concluding that RV jitter due to 
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Fig. 12. Up: Ha index vs. BIS. The dashed line is indicating the 
quiescent state. Down: H/3 index vs. BIS. For both, it is seen that 
the scatter for BIS is higher when occur flare events. Error bars 
for the indices are of order of 0.001. 



flare occurrence would be non-negligible, although probably be 
a stochastic jitter component. Chromospheric activity indicators 
exhibited an enhancement at flare state, the broad emission of 
Balmer lines and He I in emission being the most notable 
features (see Fig. [U. As pointed out by Saar (2009), although 
these lines are excluded when we measure RV, it is possible that 
a significant core filling in photospheric lines occurs when there 
is a flare event. The cause could be upper photospheric heating. 
Results by Houdebine (1992) state that heating is propagated 
down to low photospheric levels. 

A second related problem is the efl^ect of large flares on BIS. 
While it has not been studied until now, it is expected to be more 
pronounced, since bisectors are more sensitive to changes in line 
profiles. To our knowledge it is reported here for the first time. 
Fig. [12] a shows the relationship of Ha index vs. BIS, where 
the dashed line corresponds to the quiescent state, and higher 
values for Ha index indicate the occurrence of a flare event. It 
is seen that the scatter for BIS is higher when a flare occurs. 
Outliers at quiescent state correspond to a low S/N rate. Similar 
BIS behaviour is seen in Fig.[T2]b, that shows Hy6 index vs. BIS. 

5. Orbital solution for BD+20 1790 b 

We computed the orbital solution for the RV data using a stan- 
dard Keplerian fit with the RV period estimated by the Least 
squares periodogram. The fit was obtained firstly only consider- 
ing the FOCES data, averaged by night, in order to avoid intra- 
night scatter. After this, we added the SARG data to improve 
the fit. The results for the fit considering only FOCES data or 
all data from the two spectrographs were compatible within un- 
certainties. With the addition of RVs measured in winter 2008 
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Fig. 13. Radial velocity variability of BDh-20 1790,a. Up: 
Circular orbit, b. Down: Eccentric orbit. Values marked with 
circle symbol represent FOCES runs except stars that represent 
DDT FOCES 08B run. Diamond symbol are for SARG runs. 



(DDT FOCES 08b run), the Least squares periodogram is strik- 
ingly improved, and the 7.78 day peak clearly dominates the 
power spectrum. Attempts to perform a Keplerian fit using the 
second and the third highest periodogram peaks produced sig- 
nificantly worse folded curves. We have included to perform the 
fit the RV set computed by Lopez-Santiago (2005). A first fit 
(see Fig.[T3b) derives a close-in massive planet (a = 0.066 AU, 
M2 sin / - 6.54-Mjup) in a circular orbit (e = 0.05) with a rota- 
tional period of 7.7834 days and a reduced of 1.07. Also we 
present a second fit (see Fig. \T3h ) with the same period for an 
eccentric orbit (a = 0.066 AU, Masin/ = 6.15M^„p, e - 0.14, 
X'^ = 0.997). Due to the sampling of the data, we cannot discard 
a possible eccentric orbit. 

Orbital elements for both solutions are compiled in Table |4] 
and discussed in the next section. 

As additional test, we computed the orbital solution remov- 
ing the data affected by flare events. The fit derives a solution 
(fl = 0.066 AU, M2sin/ = 6.54Mjuj„ e = 0.01, K = 0.91 km 
s ') compatible with the solution when considering all the data. 



The fit is presented in Fig. [14] 
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Fig. 14. Radial velocity variation of BD+20 1790 computed con- 
sidering only the data that are not affected by flares. Circle 
symbol represent FOCES runs except stars that represent DDT 
FOCES 08B run. Diamond symbol ai-e for SARG runs. 



Table 4. Orbital Parameters of BD-i-20 1790 b 



Parameter 


Solution 1 


Solution 2 






7.7834 ± 0.0004 


7.7834 ± 0.0004 


days 


conj 

a 


3085.8 + 0.5 


3086.30 ±0.18 


HJD 


0.066 ± 0.001 


0.066 ± 0.002 


AU 


e 


0.05 ± 0.02 


0.14 ±0.04 




K 


0.93 ± 0.03 


0.84 ± 0.06 


km s"' 


J 


8.22 ± 0.01 


8.12 ±0.04 


km s"' 


(0 


200.4 ±21.8 


120.7 ± 14.0 


degrees 


M2sini 


6.54 ± 0.57 


6.15 ±0.59 




nils 


138.9 


132.3 


m s"' 


X" 


1.071 


0.997 





" Time of periastron passage 



6. Discussion 

The lack of a relation between the BIS and spectroscopic indices 
with the RV period, as well as the different RV and photometri- 
cal period strongly suggest that the RV variations are due to a 
planetary companion. However, it is possible that the RV varia- 
tions are actually due to a combination of phenomena (activity 
and planet). 

Stellar magnetic activity may be influenced and enhanced by the 
presence of a close-in giant planet, as proposed by Cuntz et al. 
(2000), Cuntz & Shkolnik (2002) and Lanza (2008). Thus, the 
presence of a planetary companion could be an interpretation 
for the high level of stellar activity detected. In a recent paper, 
Lanza (2009) proposes a new model that predicts the formation 
of prominence-like structures in very highly active stars with 
close-in giant planets. Also, as presented in Sec. 2 and Sec. 4.6, 
the large flares, with energy releases in the superflare regime, 
and the high rate of flare ocuiTence, could find a source in addi- 
tion to stellar activity in the reconnection of the stellar coronal 
field as the planet moves inside the Alfven radius of the star (Ip 
et al. 2004). In a forthcoming paper we explore in detail these 
possible star-planet interactions (Hernan-Obispo et al. 2009c, in 
prep). 

In addition, as suggested by the statistical analysis by 
Kashyap et al.(2008), the X-ray flux from stars with close-in gi- 



ant planets is on average 4 times greater than those with more 
distant planetary companions. For the 'close-in' sub-sample, the 
X-ray luminosity is Lx - 10^*^'^ erg on average. The X-ray 
luminosity of BD+20 1790 is 5 times brighter than this average 
which is consistent with chromospheric and X-ray emission in- 
duced by the presence of a massive close-in companion (Lanza 
2009). 

Even though the stellar activity could swallow the RV signal of 
a planetary companion, we can detect it for BD+20 1790 b since 
it is a massive planet. The RV variation is large enough even 
though the RV accuracy is tipically about 150 m s"' . 

Due to the observational strategy (the data are not part of a 
planet-search program), the eccentricity is poorly constrained. 
Indeed there is no "a priori" reason to discard an eccentric or- 
bit since the circularization time-scale computed is up to several 
Gyr, but more data is required to properly characterize the ec- 
centricity. RV optical and infrared follow-up over twice the RV 
period will enable us to constrain the orbital solution as well 
as confirm the presence of the planet. More massive exoplanets 
M2 sin / ~5Mjiip with orbital periods longer than about 6 days 
have eccentricities significantly larger than lower mass planets 
(Udry & Santos 2007). Another possibility is that additional un- 
detected longer period planets are maintaining the eccentricity 
of BD+20 1790 b. Both situations have been discussed in detail 
by Wu & Murray (2003). 

It is worth noting however that the star is metal-rich, as presented 
in Sec. 2. The existence of a correlation between stellar metal- 
licity and planet mass has been reported by e.g. Santos et al. 
(2001), Fischer & Valenti (2005), Guiflot et al. (2006). Massive 
planets tend to form around metal-rich stars, i.e., planets that or- 
bit around metal-rich stars also have higher mass cores. 
Compared to other planets of similar masses and orbitfl and 
taking into account statistical results described in recent reviews 
(Udry & Santos 2007), BD+20 1790 b does not exhibit unusual 
characteristics, except for its young age and its relatively high 
mass. We used a complimentary method to determine the stel- 
lar age from Mamajek & Hillenbrand (2008) (hereafter MH08), 
that uses the fractional X-ray luminosity, Rx = Lx/Lt,oi- MH08 
demonstrate that Rx has the same age-inferring capability as the 
chromospheric index By using their equation [A3] we es- 
timated an age for BD+20 1790 of up to 35 Myr. Considering a 
value for logR'^jj^ - -4.2 on average, we can also estimate the age 
with the new relation proposed by MH08, by equation [3]. We 
computed an age of up to 58 Myr. These values are in agreement 
with the range estimated by Lopez-Santiago et al. (2006). 

Lowrance et al. (2005) included this star BD+20 1790 in 
a coronographic survey for substellar companions using the 
coronograph on NICMOS/HST and the 200 inch Hale Telescope 
(Palomar Osbservatory). No companions were found beyond 10 
AU. However, the orbital solutions we find suggest a semi-major 
axis below 0.1 AU, clearly beyond their resolving capabilities. 

Great care must therefore be taken when extrapolating prop- 
erties of early stellar evolution stages from the characteristics 
of the latter stages, since the cuiTent knowledge about planetary 
system evolution is still somewhat speculative. The exoplane- 
tary zoo is such that new planets with unusual properties require 
a replanting of planet formation and migration scenarios. Planets 
discovered around young stars could be the missing link that re- 
construct the scenarios between exoplanets and protoplanetary 
disks. Indeed, further study of BD+20 1790 b has the potential 



* Observational data for the more than 370 exoplanets are compiled 
on the Extrasolar Planets Encyclopaedia \http://exoplanet.eu ), man- 
tained by J. Scheneider 
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to improve our understanding of planetary systems at early evo- 
lutionary stages. 

7. Conclusions 

This paper describes the investigation of RV variations for the 
young and active K5V star BD+20 1790. Based upon the anal- 
ysis of the BIS of the CCF, as well as activity indicators and 
photometry, the presence of a planetary companion is shown to 
be the best interpretation. The orbital solution results in a com- 
panion with a mass in the planetary regime. No photometric pe- 
riod larger than 2.8 days strongly supports the planetary origin of 
the observed RV variations. Two solutions for the orbit are com- 
puted and discussed. The presence of a close-in massive planet 
could also be an explanation for the high level of stellar activity. 
Since the RV data are not part of a planet search program, we 
can consider our results as serendipitous evidence of a planetary 
companion. Indeed additional RV optical and infrared follow-up 
will enable us to constrain the orbital solution as well as confirm 
the presence of the planet. This is thus far the youngest main 
sequence star for which a planetary candidate has been reported. 
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Appendix A: Periodogram 

We use a Least squares periodogram approach to identify and 
visually illustrate the relevant periods in the data. It differs from 
the more classic Lomb-Scargle periodogram (Scargle 1982) in a 
few key aspects. For a given period P, a linear model of the form 
Vr = y+A cos 2n/Pt+B cosln/Pt is fitted using a weighted least 
squares to the data and the of the residuals is obtained. The 
minima reveal the candidate signals of interest. One can rep- 
resent the root mean square of the residuals (RMS) with respect 
to the period to show the relevant periods as minima (This ap- 
proach is used to illustrate the photometric periods in Fig. 13). To 
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recover a more familiar view of a periodogram, one can compute 
the Power of each period P as 



that follows a Fisher-F Distribution with 2 and n„i,s - 3 de- 
grees of freedom and can be used to obtain a first hint of the 
False Alarm probability of a given solution. This definition of 
the power measures how much the^^ of the fit improves when a 
sinusoid of period P is included (see Gumming 2004 for a more 
detailed description). 

Since analytical approaches tend to give optimistic confi- 
dence levels, it is desirable to obtain the False Alarm Probability 
of a solution empirically. To make this, we generate a large num- 
ber of synthetic datasets (10^) with the same sampling cadence 
(same dates) but only containing random noise. Then, for each 
realization, we compute the Least Squares periodogram and find 
the period with higher power, which will be a false alarm. A 
histogram of False Alarms as a function of the Power is ob- 
tained and its complementary cumulative distribution gives the 
False Alarm Probability of a given peak in our signal. The Least 
Squares periodogram of the RVS data and its associated empiri- 
cal FAP probability distribution are shown in Fig.H) Compared 
to the Lomb-Scargle periodogram, this approach allows a proper 
weighting of each observation and can be easily generalized to 
include other time-dependent effects in the signal at the period 
search level. 



Power(P) 



xl/inob, - 3) 
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